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Abstract 

We consider late-time cosmology in a (phantom) scalar-tensor the- 
■ ory with an exponential potential, as a dark energy model with equa- 

tion of state parameter close to -1 (a bit above or below this value). 
k>H ; Scalar (and also other kinds of) matter can be easily taken into ac- 

$_j ■ count. An exact spatially-flat FRW cosmology is constructed for such 

theory, which admits (eternal or transient) acceleration phases for the 
current universe, in correspondence with observational results. Some 
remarks on the possible origin of the phantom, starting from a more 
fundamental theory, are also made. It is shown that quantum gravity 
effects may prevent (or, at least, delay or soften) the cosmic doomsday 
catastrophe associated with the phantom, i.e. the otherwise unavoid- 
able finite-time future singularity (Big Rip). A novel dark energy 
model (higher-derivative scalar-tensor theory) is introduced and it is 
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shown to admit an effective phantom/quintessence description with a 
transient acceleration phase. In this case, gravity favors that an ini- 
tially insignificant portion of dark energy becomes dominant over the 
standard matter /radiation components in the evolution process. 
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1 Introduction 



Recent astrophysical data, ranging from high redshift surveys of supernovae 
to WMAP observations, indicate that about 70 percent of the total energy 
of our universe is to be attributed to a weird cosmic fluid with large and 
negative pressure, the dark energy (see [1, 2] for a recent review) and that 
the universe is currently in an accelerating phase. It also turns out that 
the dark energy equation of state parameter w is close to —1. So far, the 
simplest possibility proposed for such kind of dark energy is the use of a 
scalar field (or a scalar-tensor theory). However, scalar-tensor theories are 
not free from problems, especially when they are considered directly as dark 
energy candidates. 

Much attention has been drawn by scalar fields in studies of the early 
time universe. A variety of scalar potentials have been considered and a 
number of accelerating (inflationary) cosmologies have been advocated. For 
instance, the interesting quintessence model [3] with w slightly bigger than 
— 1 is quite popular for the explanation of early (and late) time acceleration, 
especially in the case of exponential potentials [4]. Moreover, exponential 
scalar potentials often appear naturally after compactification in string/M- 
theory. Needless to say, such description is model-dependent and is still quite 
far from the final goal: the formulation of a plausible and consistent dark 
energy theory. 

Another line of research is related with the case where the dark energy 
equation of state parameter is less than —1, since this possibility is not ex- 
cluded by astrophysical data. The typical example of a dark energy of this 
kind is provided by a scalar field with negative kinetic energy, dubbed phan- 
tom (see [5, 6] and references therein). At first sight, such models may look 
rather strange and they lead to a number of unpleasant consequences, as a 
finite-time future singularity (the Big Rip) [7, 8, 9]. Nevertheless, the pos- 
sibility of negative energies seems to be acceptable in classical scalar-tensor 
theories. Actually, many of them do contain phantoms, as the ones coming 
from string/M-theory compactification, or higher-derivative super gravities, 
or modifications of Einstein gravity itself. In fact, the issue is somehow del- 
icate, since what looks like a phantom in one reference frame may radically 
change its nature in another frame (e.g. after a conformal transformation). 
In this sense, even in the absence of fundamental, physical meaning the phan- 
tom can be still useful as a convenient mathematical tool in order to study 
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cosmological models in standard scalar-tensor theories, because a phantom- 
related frame may lead to a simpler formulation of the problem. Finally, 
there are examples where an effective phantom/quintessence description of 
the late-time universe naturally appears, even if the starting theory does not 
explicitly exhibit the phantom/quintessence structure. 

In the present work we study different cases of a late-time spatially- flat 
FRW cosmology in the (phantom) scalar-tensor theory, mainly with exponen- 
tial potentials. Such scalar is considered as a dark energy and the possibility 
of deriving the current speed-up is shown also in the presence of matter. Ex- 
act FRW cosmologies are constructed for the (phantom) scalar-tensor theory 
with an exponential potential, a model that can be important for under- 
standing attractors and the stability properties. The possibility of avoiding 
the unwanted Big Rip by simply taking into account quantum gravity ef- 
fects, which may become dominant near future singularity, is demonstrated. 
Finally, the present dark energy dominance and acceleration, within the ef- 
fective phantom/quintessence description, is discussed in the model where a 
new form of higher- derivative, gravity-matter coupling is introduced. 

The organization of the paper is the following. In the next section 
spatially-flat FRW cosmological solutions are discussed for scalar-tensor grav- 
ity with scalar matter. Explicit examples of accelerating (and decelerating) 
scale factors are presented for exponential potentials when the theory con- 
tains one or two scalar fields. In Sect. 3 the general solution for a spatially- 
flat FRW cosmology which includes eternal or transient acceleration is found 
in the (phantom) scalar-tensor theory with exponential potential. This is 
based on the use and extension of a method recently developed by J. Russo. 
The comparison with particular solutions of the previous section is done. 
Sect. 4 is devoted to the study of the influence of quantum gravity effects 
on the Big Rip singularity in phantom cosmology. It is shown that taking 
them into account properly may change the future of the universe, from that 
with a finite-time singularity to an ordinary deSitter space. In Sect. 5 a 
higher derivative matter-gravity coupling is suggested as a new sort of dark 
energy model. It admits an effectively phantom/quintessence description 
and does explain the current dark energy dominance over standard matter 
by gravity assistance. Stability analysis of the model demonstrates that the 
acceleration phase is actually transient. A summary and outlook are given 
in the discussion. In the Appendix we outline how the (phantom) scalar- 
tensor theory may originate, via compactification, from a higher-dimensional 
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(super)gravity. 



2 Examples of the accelerating universe in 
(phantom) scalar-tensor theories. 

We start from the action of multi-scalar-tensor theory. Several illustrative 
examples of the (accelerating) FRW cosmology will be presented here as 
simple dark energy models (see [1, 2] for a recent review). A scalar field, </>, 
which may be later regarded as a phantom, couples with gravity. As is typical 
in these models, a second scalar field, x, is considered. The string-inspired 
Lagrangian in the d-dimensional spacetime is: 



Here a and 7 are constant parameters and V((f>) (U(x)) is the potential for 
4> (x)- If the constant parameter 7 is negative, <fi has a negative kinetic 
energy and can be regarded as a phantom [5, 6]. We should note, however, 
that 7 need not be negative in order to obtain the accelerated universe as 
we will later see. As the matter scalar x does not couple with directly, 
the equivalence principle is not violated, although the effective gravitational 
coupling depends on as . One may go to the Einstein frame by the 

scale transformation 



In the following, the quantities in the Einstein frame are denoted by the 
index E. After the transformation (2), the action (1) has the following form: 




(1) 



9ixv = e d ~ 2V g E ^u ■ 



(2) 




(3) 
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In the Einstein frame, there appears a term coupling the matter \ with (p. 
Even if 7 is negative, when 

<^ + i>0, (4, 

the kinetic energy of becomes positive, as for a usual scalar field. Hence, 
the remarkable observation follows that, what is a phantom in one frame may 
not be a phantom in a different frame, especially if the coupling is taken into 
account. 

As first step, one considers the d — 4 case and assumes x — 0- We now 
define (p and V(ip) as 



<p = ^a* + l, V( ( p)=e- a *V(<l>) , (5) 
and assume that the metric has the FRW form in flat space: 

ds 2 E = -dt 2 E + a E (t E f ■ ( 6 ) 

i=l,2,3 

Here t E is the time coordinate in the Einstein frame. When ip only depends 
on the time coordinate, the FRW equation and ip equation follow: 

Here the Hubble parameter (in the Einstein frame) is defined by H E = ^^j- 
Note that models of such sort may have a double interpretation: as multi- 
scalar-tensor theories or as matter-scalar-tensor theories. In other words, 
some scalars may be considered as matter or as part of a gravitational theory. 
It has been suggested that such models may describe the inflationary early 
universe as quintessence [3]. 

Special attention in cosmology has been paid to exponential potentials [4], 
which often follow from string/M-theory compactification. If V(<p) behaves 
as an exponential function of ip 

V(tp) ~ V e~ 2 ^ (9) 
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(Vo and (p are constants) during some period, as in the present universe, the 
solution of (7) and (8) exists: 4 

a E = a E0 [ — ) , (p = (p m—. (10) 

Here 



In the original (Jordan) frame (1), the time coordinate t and scale factor a(t) 
in the FRW form of the metric 

ds 2 = -dt 2 + a (t) 2 ( dxl ) 2 > ( 12 ) 

i=l,2,3 

are related with the corresponding quantities in the Einstein frame (6) by 
dt = e-i*dt E = ^f^d (13) 



4 The action (1) with the potential (9) belongs to the class discussed in [10]. In fact if 
we redefine the field cj> by 

the action (1) can be rewritten as 

+ 1 <fV^(-^U^x-t%)) • 



F($) = — <f 2 , I7(*) = Vb ' 1 



Here 



4 n ' v ' v 2 vh; 

Then, one can use the same arguments as in [10] in order to fit the parameters so that they 
satisfy the present cosmological data. The non-minimal scalar-gravitational coupling term 
which is required by renormalizability of the quantum field theory in curved spacctimc 
[11, 12] may have very interesting effects on the phantom cosmology (see [13] for a recent 
discussion.) 
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a<f> 



a = e 2 a E = a E0 [- — I = a E o I — I (14) 



,3,oV ln r- (15) 



Here 



/^ — =g=, t = ^-. (16) 



Then, cosmic acceleration (d > and a > 0) occurs when 

£ 



3, 2 _ <A)Q 

, >!■ (!7) 



By properly choosing the parameters a, ipo, and 7, the present cosmic accel- 
eration can be realized. For the matter with p = wp, where p is the pressure 
and p is the energy density, 

a oc tsp^Fij . (18) 

Then, effectively 

_ 1 1 3 V 2 v / ^+iy / ¥?o (2/3-9^+8 
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As w diverges when the denominator in the second term vanishes, w can take 
any value by properly choosing a, ip , and 7. For example, if 



(p a 



w — — 1. Several interesting cases deserve attention: 
when (fl > - 

w<-i if |^>-»r>2 



(20) 



2 ™ v^+f 

yo« v. 3,.2 ~_ <gog 



?/ , > _i if yo Q > 3 2 g < 9 



(21) 
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when < - 

w<-l if 2 > -g^= > 
w>-l if -^=<|^or-^=>2. 

As is clear from (18), cosmic acceleration (a > 0) occurs when 

u> < - J . (23) 

When w < — 1, the universe is not expanding but shrinking, although the 
universe is still accelerating. In the case w < — 1, if we change the direction 
of time as t — > £ s — £, the universe is both accelerating and expanding. 
In the especial case when w = — |, we find 

^ = | ■ (24) 

Note that even when 7 is positive — e.g. the kinetic energy of </> is positive 
as for usual matter — the effective w can still be less than —1. For example, 
with the choice 

^0 = 4 , ^ = 3 , (25) 
it follows that (3 — and from (19), 

UV2 + 203 

213 v ; 

Note that the above considerations are appliable both for the early as well 
as for the late time universe. Indeed, the simple solution under discussion 
(especially, the no phantom case) has been considered in different situations. 
Later on, the general solution for the scalar-tensor theory with an exponential 
potential will be discussed. With a known general solution, it is much easier 
to understand the type of situation which appears: if it is a transient (or 
eternal) acceleration, or an attractor, or something else. 

In the previous example the matter field \ is zero (x = 0). We now 
consider the case of x which is much related with the so-called double 
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quintessence model [14]. In the Einstein frame (3), the FRW equation, <p 
equation, and \ equation have the following form: 



-2(3K 2 e-^U{ X ) , (28) 

Q = -L( e-^^-) + 3H E e~^^ + e~ 2 ^U\x) • (29) 
Me V J Me 

We again assume the exponential potentials 

V(<p) = V e~^ , U( X ) = U oX 4 -^ ■ (30) 
The following Ansatz exists 

H E = , V = Vo In r^- > X = Xo(-M • (31) 

The parameters /io, £eo, and xo should be a solution of the following algebraic 
equations, which can be obtained from (27), (28), (29): 

= -S/ig + ^ + ^ + ^^ + y^xr' 4 , 
= 3(-^ + 3Wo)-^ + ^T^-2^o4oXo"^, 
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= _ T ^ T + 1 j + ^ +4[/A ^__^ o »„. (32) 

For instance, for the special example 

[27 

V?o = ^/§, / 9 = 2 y I 7 ) ' y o = ' ( 33 ) 
the explicit solution follows: 




5 V5 HT 



Using (16) and (33), one arrives at 

4Q 

7 = -^« 2 . (35) 

Then is surely a phantom with negative kinetic energy in the physical 
Jordan frame (1). Since U(x) = UoX 2 > U(x) corresponds to a mass term and 
the mass m x is given by 

m\ = 2U . (36) 

Also 

'» = S£ ■ (37) 

Since dt = ±e~^dt E = ±e~^dt E = ±^-dt E , we find 

— = ±ln^. (38) 
Then, in the physical Jordan frame one gets 
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a = e 2 a E = a E o [ - — = a E0 e 12t Eo , (39) 



_a± ( t E \ 

, P 2 , t E , 2 t 

4> = -Lp = -In— = ± , (40) 

a a t E Q cet E o 

X = xJr~ = Xoe ±T ^ . (41) 

Since the scale factor a behaves as an exponential function of t, the Hubble 
parameter H is a constant 

H = T ^ = T 7 -^, (42) 

which is linear in the mass m x of %. Eq. (42) also tells that the effective w 
is —1 (similar to the cosmological constant). 

As a different, more complicated example one can consider the case when 
a = 0, and therefore (3 — 0, but the potential depends on both <fi and %: 

S = J d A x^~g [^(R- \d^^j - \d»xd»X ~ W(<p, X ) } • (43) 



10 



As a = 0, the Einstein frame can be regarded as a truly physical one. Then 
the equations corresponding to (27), (28), and (29) are: 

** - i© 2 + ?(§) 2 + ^^ (44) 

° = 3 (§ +3 40 + ^*)- <«> 

Here the derivative of W(ip, x) with respect to ip (x) is expressed as W jtp (ip, x) 
(W jX (</?, x))- The case of the exponential potential W(<p, x) ma Y De of interest 

W(<p, x) = W e^ x "~^ , (47) 
where Wo, rj, xo, and ip are constant parameters. Assuming 

h t t 

H = -j- , up = ip In — , x = Xo In — > ( 48 ) 
i r to 

with constants h and t , Eqs. (44), (45), and (46) reduce to the following 
algebraic equations: 

fl 4. ^ 4. ^ 4- 

U - 6n o + -^ + — + —2 — 

= (-1 + 3/^) -(2 + 77) K 2 Wb*g, 

= 3 X l(-l + 3h ) + W 2 W tl . (49) 

Eqs. (49) give 

- 2 xl = , (50) 

ho = (51) 
2(2 + 77) 

9 W ,2 3 ^0 ( 3/ *0 - 1) /r X 

W = • ( 52 ) 

-2 < T) < . (53) 



Eq. (50) shows that 
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The effective w is given by 



2 4(2 + 77) 

W = -l + — = -1+ l Q / 

3/i 9y9g 



(54) 



As discussed in (53), as w > — f we do not have a phantom. If the second 
term in (54) is small, one may obtain quintessence. 

Other types of matter can be easily considered too. For instance, matter 
may be dust. The energy density pdust in the Jordan frame (1) behaves as 
Pdust oc a~ 3 . Then 



Pdust = Poa = Poe 2 = p e 2 ^ , 



Pdust 



= . 



(55) 



It is well known that dust has no pressure. It could correspond to the baryon 
and/or cold dark matter components. Instead of (27) and (28), we obtain 
the following equations of motion in the Einstein frame: 



3 f dip 

4 \dl~E 



1 ~. 
2 

dip 



K 

~2 



dt'g dt e 



(5k 2 



+ V\ip)-^-e~^p af 



With the form of V(<p) as in (30), the solution occurs 



, ^ f tE 

ip = <p m — , a £ = a — 



2_ £-£0 
3 6 



Here 



r - 



/5Vo + 9yg - 9/3^o 



6^o 

Po K 2 e-^(/?Vo + 9^-9^o) 



2F (8 - 2/?<p - 9<pg) 
The time-coordinate t in the Jordan frame is given by 



t 

T 



tE_ 

tEO 



to 



(56) 
(57) 



(58) 



(59) 



(60) 
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and the scale factor can be obtained: 



2 2(1*0) 
/ I \ 3 2-/3^ 

a=ao UJ ■ (61) 

The effective w is found to be 

W = ~ l + W^)- (62) 

Hence, if 

1 < (3(fo < 2 , (63) 

k; < — 1. If we assume V > (and £q > and po a o 3 > 0)> Eqs. (59) require 

/?Vo + 9c/? 2 - 9pVo > , 8 - 2p> - > , (64) 

which give /3 2 y?o — UpVo + 8 >, that is, 

_ ll->/89 _ ll + v 7 ^ 

pVo < g or > 2 ' ' ' 

Numerically, pVo < 0.7830... or f3ip Q > 14.9339..., which contradict the results 
in (63), and there is no accelerating universe. With the assumption V < 0, 
the accelerating universe takes over. For instance, with the choice fiipa = § 
and (pQ = | (P = ±| and po = ±|), we obtain t 2 , = — 2^ and a[j = 

P0« C 2 



2V 

Going back to the case of two scalars, by variation over g^ v the Einstein 
equation follows: 

it = ^ (-\d P <t>d p <j>g» v + id^d v <\> - V(4>) 9lH , 
+2e- a<t 'V lx V u (e^) - 2<^e"^V 2 (e^)) , 

T* u = e~°* [~\d pX d P X9^ + d»xd„X - U( X )gJ) ■ (66) 
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T^ u and T* v may be regarded as the effective energy-momentum tensor of 
and respectively 5 . In particular, in the FRW metric (12), we find 

Tt = P 4, 

= ^{^ 2 + V(d>)-6aH^ , 

n = 7 



Tf. = p*a% 



V((f)) + 2a<p + 2a 2 <p 2 + 4aH<p} a 2 ^ , 



The effective «r and w x can be defined as follows: 



p- = 



. 2 ' 



(67) 



w+ = 4 , w x = — ■ (68) 
p<t> p x v ' 

For the first example in (13), (14), and (15), it follows that 
3 (2/3 - 3y9 ) \l 



* W ~ 3^o) {Vl gg - 9y9 ) + 8^0} 
8tff*(l-%a) 2 

which agrees with w = - = — in (19). Note that the consideration of 
various entropies for dark energy models can be done and then interesting 
holographic relations among them occur (see [15] for a recent discussion). 
For the second example in (39), one gets 

4> - 14 4> - 19 



5 The usual energy-momentum tensors Tfy x given by the variation over g^ v are related 
with T/fc* by T$* = e a <t>T$>x . 
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Hence, 

™* = ~, ™ x = \- (71) 

However, since 

147 

P* + P X = -(P* + P X ) = ^7t- , (72) 

it turns out that to = — 1, which is consistent with (42). We should also note 
that, in the Einstein frame (3), (j) has a positive kinetic energy. 
For the case of (48) with (50), (51), and (52) we find 

p= ^ + ^- w(M = &^»%, (73) 

which reproduces (54). Having these various examples of the (accelerated) 
evolution of the current universe one can compare it with recent astrophysical 
data in the way discussed recently in [2, 16]. Of course, above illustrative 
examples of current speed-up may correspond to transient acceleration. In 
other words, stability of the solutions pretending to be realistic ones should 
be investigated in detail. 



3 Exact FRW cosmology for the (phantom) 
scalar-tensor theory with an exponential 
potential. 

In the present section the exact FRW cosmology in the (phantom) scalar- 
tensor theory with an exponential potential will be discussed. First, the 
method which is appropriate to obtain the exact FRW solutions will be re- 
viewed [18] (for introduction of similar variables in quantum cosmology, see 
[19]). Subsequently, the formulation is extended to the phantom case with 
negative kinetic term. 

The action of a scalar field coupled with gravity is: 

S=^J d'x^-g (R - \d^d^ - V(</>j) . (74) 



15 



This action can be regarded as the action with a = in (1) or that obtained 
by replacing + \ an d e _3=5, V(0) in (3) with | and V ((/)), respectively. 

For the standard scalar 7 > and one can normalize to be 7 = 1 but for 
the phantom field with negative kinetic term, we have 7 < 0. 
For the FRW metric 

ds 2 = -dt 2 + a(t) 2 ( dxi ) 2 > ( 75 ) 

i=l,2,3 

the action (74) can be rewritten as 

S = -^ J d 3 xdt |-6ad 2 + a 3 (|</> 2 - V(<j>)} } . (76) 
The potential V{4>) is chosen to be 

2<f> 

V{(f>) = V e~o , (77) 

with constants Vb and 0o- 

First we review the standard case with 7 > following to [18]. The field 
variables a and are written in terms of new fields v and u as 



and a new time variable r is defined by 



a = e 3 , = v ' , (78) 



dr = dt\ /— ^e . (79) 



Then, the action (76) acquires the following form: 

dv du 



e „ +u _ s( „_ u) ^ _ ^ _ ^ 

0OA/37 



Varying over i> and it, the equations of motion follow: 

d 2 u , . / du\ 2 , . 

= ^ + < 1+5 >U) - (1_a) ' (81) 

d 2 v _ N / du^ 2 



= _ + (!-«)»_) _(l + a) . (82) 
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Since the Hamiltonian H conjugate to r is given by 

I fW 



H= — 





dv du 


1 d 3 x 




_dr dr 



v+u—a(v~u) 



(83) 



the Hamiltonian constraint H = yields 

dv du 
dr dr 

In terms of the new variables V and U, which are given by 

V = e (1 "" ),; , U = e (1+s)u , 



the equations of motion are 
d 2 U 



d 2 V 



When \a\ < 1, the solution of (86) is given by 



U = u + c 



, V = v+e 



-Vl-a' 2 _|_ v g-rVl-a 2 



(84) 



(85) 



(86) 



(87) 



with constants of the integration u± and v±. The Hamiltonian constraint 
(84) restricts the constants as 



u + v^ = —u_v + . 
Then the spacetime metric has the following form [18]: 



(88) 



ds 1 



3V 



+ V-e 



dr 1 



a 2 ^ 3(i-a) 



(89) 



i=l,2,3 
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When r -> +00, Eqs. (78) and (79) give 



a/1 _ ^2 / o a a a 2 r 

t - w+^v^v^^^ 2 ^ 67 ^- (90) 

2 

Here t + is a constant of integration. Hence, a oc t^ 1 and the universe is 
accelerating > if 

a 2 < i . (91) 

Note that if t>_ = n_ = 0, the behavior in (90) is exact even if r is not large. 
The case with i>_ = w_ = corresponds to the solution in Sect. 2. On the 
other hand, when r — > — 00 

1 1 _ 2t 

3(i-a) 3(i+a) ., A — tt 



* - *o- - J^I^I^e"^ . ( 92 ) 

cr V -jI/o 

2 

Here to- is again a constant of integration. Thus, we find a oc (— then 

the universe is shrinking but accelerating > 0^ if (91) is satisfied. To 

summarize the \a\ < 1 case, if V- = U- = we find an eternal expanding 
solution as in Section 2. In the general case, the solution is a bouncing 
universe, where first the universe shrinks and, after that, it expands. 
When|a| > 1, the solution of (86) can be written as 

U = u c cos (t\/ a 2 — 1 j + u s sin (t\/ a 2 — 1 j , 

V = v c cos (t^/ a 2 — 1 j + v s sin (t^/ a 2 — ij , (93) 

with constants u c , u s , v c , and v s , which satisfy 

v c u c + v s c s = . (94) 

The solution of (94) can be given by means of three independent parameters 
v , u , 9 as 

v s = v sin9 , v c = Vqcos^o , u s = UqcosOq , u c = —u sin6 . (95) 
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As a result, Eq. (93) simplifies 

V = v cos (Wa 2 - 1 - O ) , U = u sin (V Va 2 - 1 - O ) • (96) 

As 6*o can be absorbed into the constant shift of r, in the following we choose 
#o = 0. The spacetime metric looks as: 



8 2a 2a 



cos !-« ^rV a 2 — 1 j sin^" ^rV a 2 — 1 j dr 2 
+v^ E1 ufi cos^T (Wa 2 - l) sin^y (r^/a 2 - l) . (97) 



There are singularities when 

rV a 2 — 1 — nn , or + 7r . (98) 

Here n is an integer. If we write r as rV a 2 — 1 = n7r + 5r and assume 5r is 
small, we find, by neglecting numerical factors, 

t~(5r)^ , a ~ (5t)3TtW ,v, t3(2< i 1 + i) _ (gg) 

Note that > as | or] > 1. Then r = corresponds to t — 0. At t = 0, 
the size of the universe diverges (vanishes) when 2a + 1 < (2a + 1 > 0). 
On the other hand, if we write r as T\Ja 2 — 1 = (n + |) 7r + 5r and assume 
5r is small, we find, again neglecting numerical factors, 

t ~ (5t)^ , a ~ (^r) 3 ^) ~ t3 ( i-2«) _ (100) 

Note that 1 ~ 2 f_ 1 ' 1 > once more, and r = corresponds to t — 0. Then, at 
t — the size of the universe diverges (vanishes) when 1— 2a < (1— 2ct > 0). 
These results for the standard scalar with 7 > are given in [18]. 

We now extend the above formulation to the case of a phantom with 
7 < 0. For this situation, we define a complex field z and its complex 
conjugate z* by 

a = e 3 , </>= , (101) 

and define a new time variable r as in (79) by 



/3Vn 2i(z-z*) 
dr = ±dt\/-^-e . (102) 
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The action (76) becomes: 



K 



Here 



dz dz* 
dr dr 

4 



+ 1 



1 z+z*—ia(z—z*) 



a = 



(103) 



(104) 



The sign =p in (103) corresponds to the sign in (102). Varying over z*, one 
obtains the following equation: 



d 2 z , ~ x ( dz \ 

o = S3 + (i-«) ( Tt ) -d- 



Now, the Hamiltonian H conjugate to r is given by 



K 



1 8V 





'dzdz* ^ 


1 d 3 x 




dr dr 



^v+u—a(v—u) 



and the Hamiltonian constraint has the following form: 

dz dz* 



= 1 . 



dr dr 

By defining a new variable Z as 

Z = e (1 "^ } 

Eq. (105) can be rewritten as 

d 2 Z 
dr 2 

The solution of (109) is given by 

Z = z + e 



= ( 1 + a ) Z . 



(105) 



(106) 



(107) 



(108) 



(109) 



(110) 



with complex constants of integration z±. The Hamiltonian constraint (107) 
restricts the constants to satisfy 



Z + Z_ = —Z-Z^ 
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111) 



By using three real independent parameters b± and 6>o, the solution of (111) 
is given by 

z+ = b + e te ° z_ = ib„e~ ie « . (112) 
By using z± in (112), the metric of the spacetime has the following form: 



ds 2 = 



I / ~ J 1 — ia 



dr 2 



z + e 



I ^2 I ^2 I 3(l-ia) 

"V i+a _|_ z ^e~ T ^ 1+ 



E (^) 2 • ( 113 ) 



i=l,2,3 



When r — > oo, from (101) and (102), it follows that 



a ~ e 



TF + i(I% R {( 1+ ^ ln2+ } 



t ~ t + =p — J e i+a* u ; +J e vI+S^ . (114) 



2 w 3F 



Here t+ is a constant of the integration. Hence, a oc (=F (t — t + )) ss? , which 
tells that the universe is accelerating since a oc |a 2 ^|« 2 + 1 j (± (i — t+)) 2 > 
0. The effective w is given by 

u, = -1 -a < -1 . (115) 

The case corresponds to a phantom. In general, 



Q T/l 1 _ 9 

q _ ^3(i-is) l-ia ^*3(i+iS)^i+i 

8 

2 (1 + is) (2 - 3ia) 
9(1 -m) 



4 ( 1 + a 



9 



(116) 



2 2 (1 - ia) (2 + 3ia) , , 2 
+ 9(l + ia) + " . 



If \Z\ is large, then a > 0. We should note that when Z- = 0, Eq. (114) gives 
an exact solution corresponding to those in Sect. 2. In this case, Eq. (114) 
is valid even if \t\ is not small. The solution can be regarded as an attractor. 
In fact, when t ~ t+, all the solutions behave as this one. 
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On the other hand, when r — > — oo one gets 

— / 2 -2\ H( 1 ~ i§l ) lnz -} 
a ~ e zVT ^ 3 ( 1+a ) 



t „ t_± v 1-1 - . _ e -I^"l^^--/ e v^ , (117) 



3V 



Here t- is a constant of integration again. Thus, a oc (± (t — t-)) and, 
once more, the universe is accelerating. 

The solution (113) is almost given by the analytic continuation a — > ia 
of that (89), which corresponds to the standard (non-phantom) scalar. The 
behavior obtained for the solution (113) is, however, rather different from 
the non-phantom case. In the case of a non-phantom with 7 > in (89) 
—which has been investigated in [18] — when \a\ < 0, the behavior in (90) or 
(92) shows that there is a singularity only in the infinite future or past, since 
r — > ±00 corresponds to t — > ±00. On the other hand, when \a\ > 0, the 
behavior in (99) or (100) indicates that there might be Big Rip [8, 7, 9] or Big 
Crunch singularity in the finite future. Even for non-phantom matter, when 
the strong energy condition is applied, a finite-time future singularity may 
occur[24]. In the phantom case in Eq. 114) or (117), r — > ±00 corresponds 
to t — > 0. The singularity occurs in the finite future or past. Thus, if the 
universe is expanding, there should be a finite-time future singularity. This 
singularity is nothing but the Big Rip [8, 7, 9]. If there is a singularity in 
the past, the universe is not expanding but shrinking. In this sense, the 
solution with a past singularity is related with that with a future singularity, 
by reversing the direction of time. 

It is possible to relate the action (74) with the action (3) in the Einstein 
frame by identifying g^ u , 7, and V{<p) with gE^v, ( w ^h d = 4), and 

e~^~2^V((/>) (with d — 4 again), respectively. The physical metric is obtained 
by rescaling the metric as the reverse of (2). Instead of (75), we now assume 
that the physical metric is given by 

ds 2 = e^(-dt 2 + a(t) 2 ( dxl )j » ( 118 ) 

V i=l,2,3 / 

For r — > ±00, behaves as 

^T -- * hi\t-t±\ . (119) 
«V-37 
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2 

Since a behaves as a ~ \t — t±\ ^ , if 

P=-lfl, (120) 

the singularity corresponding to r — > +00 is cancelled but there remains a 
singularity corresponding to r — > — 00. On the other hand, if 




(121) 



the singularity corresponding to r — > —00 is cancelled but there remains a 
singularity corresponding to r — > +00. In the metric (118), the cosmological 

time t is defined by dt = ±e^dt, then in case of (120), we find t oc t + 

2 +1 

1 1 — t ± 1 . Here t is a constant of integration. The limit t — > £± (r — > 
±00) corresponds to i — > t . On the other hand, in the case of (121), t oc 

+ \t — t±\~sS I+1 . Here is a constant of integration. If a 2 > |, t — > t ± 
(r — > ±00) corresponds to t — > t 0) again. If a 2 > |, however, t t± 
(r — > ±00) corresponds to t — > ±00. Hence, the singularity does not occur 
within finite time. This example shows that the the type (or even the presence 
itself) of the singularity is also related with the choice of physical metric 
(frame). 



4 Quantum effects may change the finite time 
future singularity. 

Let us again start from the scalar-tensor theory with a single scalar which 
can be an effective phantom: 

L=i^+^^-^)j , (122) 

where 7 = ±1. It would be interesting to investigate the quantum properties 
of such scalar-tensor gravity. Indeed, it is known that the phantom theory 
develops a catastrophic instability at the quantum level. Hence the point is 
that taking into account quantum gravity effects (or, simply quantum effects) 
could improve the situation. 
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The calculation of the one-loop effective action in the former, 
non-renormalizable theory may indeed be performed (using the above 
parametrization and some choice for the gauge condition). The result is 



loop 



+ 



7 



V - 2V" 



d A x, 



\v*-HV? + \(v»y 



-g 

3 12 



+ 1^2 _ l R( p 

40 6 Y " 



(123) 



The above one- loop action is found in Ref. [20]. In order to consider this 
effective action as a finite quantum correction to the classical one, the cut- 
off L should be identified with the corresponding physical quantity. For 
instance, when the universe is in the (almost) deSitter phase, the natural 
choice is L 2 — \R\, as the curvature is strong enough and constant [11, 12]. 
At the same time, in the region where \V\ 3> \R\, L 2 should be identified 
with |V|. 

Hence, even in the situation with V — 0, starting from the action (122) 
with the usual scalar, the phantom terms may be induced. This happens 
if the universe goes through a region with negative curvature. With ac- 
count to the potential, for some fine-tuning of V one can again arrive to the 
QG-induced phantom theory, which subsequently can change the universe 
evolution. 

Here, we consider the action where L 2 is replaced with \R\ as a simple 
example: 



loop 



R\ 



d xy— ^ln — — < -V 



l(V'f + -(V") 



1 



ll\2 



+ 



2V" 



3 12 



43 2 
R + 60^ 



(124) 



The variations of this action are given by 

1 ^V^l-loop 



2 /i 2 



2 



7 iV' f + - (V") 



ll\2 
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+ 



7 



V - 2V" 



- 2V M 



7 



V - 2V" 



3 12 



7, 



1 tfwi-loap _ 


2 /i 2 


. 2 I 2 


7 (Vf + ^ {V"f 






+ 


-V - 2V" 
2 




3 12 


R + m R ^ 


+ — R 2 -^R(j)a^ + - 

40 6 4 








ly - 2V" 
2 


+ 


3 12 


R^ 



- (v^V - g^ u V 2 ) 



13 
~3~ 



43 

R»P 

30 ' 



43 



+— {(V^i^ + V a V»R av ) - V 2 R pv - g^V p VoR p(T } 

+—RR^ U + — (VV - g^V 2 ) R + ^R^6 > p 
20 20 v y ' 6 ^* p ^ 



7 



+ (-IF" + VV - ^V 2 ) 



1 

2R 



+ 



7 



V - 2V" 



3 12 



;V 2 -7(^ 

„ 43 „ 
R+—R 
60 ' 



/\2 



1 



//\2 



l Q/3 



R 
(125) 



(126) 



In the case of occurrence of the Big Rip singularity, the curvature quickly 
grows. However, this means that quantum effects (e.g. quantum gravity 
effects) become important not only for the early universe but also for the 
future universe. These quantum effects may even become dominant when 
the universe approaches the Big Rip. Suppose that the quantum correction 
becomes dominant owing to the fact that Vl^i-ioop contains higher derivative 
terms. In this case one can neglect the classical terms. To simplify the 
situation even more, we assume that the curvature and the scalar field <fi are 
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constant 



D _ 12 



The potential V(d>) is chosen as the exponential function of (ft: 



(127) 



V(<f>) = He- 
Then from (125) and (126), we obtain 

1 SWl-Uyap 



= 



f=g 5(ft 
1, \R\ 

2 111 7^ 



4/5 4j 8 
"0o~ \2~ 0f + 



K e *° 



= 



2 / 13 7 

1 W-loop 



-iln 

4 



12 



5 47 
2~ 0^ + 



13 7 
T + 30| 



T , _2£l2 147 

u,e ^0 1 

Z 2 5/ 4 



+ 



H V 



(128) 



(129) 



5 47 
2~0f + 



+ - 



13 7 \ Tr _2c 441 



2/ 2 V 3 ' 34ft J ' u ~ 40/ 4 _ 
Eq. (129) can be solved with respect to I 2 : 



(130) 



'5_47 _8\ /13 1 
' 2 ^ + 0^ V 3 30g 



V e *o . 



(131) 



We should note, however, that Eq. (130) is not consistent with the expres- 
sion in (131) in general. Then, Eq. (130) might be regarded as an equation 
determining /i. We should also note that the r.h.s. in (131) is not always pos- 
itive. In the case 7 > 0, when 7 2 < 5, the r.h.s. in (131) is positive, but when 

4 



7 2 > 5, it is positive if 



7 + V7 2 - 5J or (ft 2 <l[l- V7 2 - 5 j • 0n 
the other hand, in a phantom case 7 < 0, the r.h.s. in (131) is positive if 
0o > — yjt Anyway, there may occur a (asymptotically) de Sitter solution. 
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Thus, before entering the Big Rip singularity, the universe becomes a quan- 
tum de Sitter space. This qualitative discussion indicates that the finite time 
future singularity may never occur (or, at least may become milder) under 
the conjecture that quantum effects become dominant just before the Big 
Rip. Due to the sharp increase of the curvature invariants near the Big Rip, 
such a conjecture looks quite natural. 

A similar phenomenon occurs even without quantum gravity. Indeed, let 
us consider again the phantom theory of Sect. 3 with the same potential. 
For the FRW background, if we assume that 4> only depends on time, the 
equation of motion for <j) is given by 







The energy density is 



"7 



P<t> 



7 
2 



dt 



+ V{<j>) , 



and the FRW equation has the following form: 

6 



:H 2 



K 



Here H = Then a solution of (132) and (134) is 



t 9 -t 



tl 



H 



7« 



4 (t s - 1) 



= 0o In 

Here t s is a constant of integration and t\ is given by 

i4>l ( i - ¥ 



t 2 



Eq.(135) shows 



2Va 



t,~t 



7« 
4 



(132) 



(133) 



(134) 



(135) 



For a phantom with 7 < 0, a grows up to infinity at t 
Rip singularity [8, 7, 9]. 



(136) 
(137) 

t S) which is the Big 
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In general, Eq.(132) shows that 



which is positive if 7 < 0, i7 > 0, and <j> 7^ 0. Hence, for the phantom with 
negative 7, the energy density increases in general. We should also note that 
since the contribution to p^ from the kinetic term is negative if 7 < 0, we 
find 

P4> < V{<t>) . (139) 

For the case V(<f>) = 0, the energy density p^ is not positive. In general, the 
Big Rip singularity occurs due to the rapid increase of the energy density of 
the phantom scalar. When V((f>) = 0, the singularity does not occur. To be 
concrete, we also consider here matter to be dust, whose energy density is 
given by 

Pd = g , (140) 

with a constant p , which we assume to be positive. When V((f>) = 0, the 
solution of (132) is given by 

Here c is a constant. Then the FRW equation has the form: 
This equation (142) can be solved easily as 



a 3 = -^— + -^{t-t s ) . (143) 



2p 4p, 



Here t s is a constant of the integration. Then, there is no singularity and 
there is not acceleration, either. 

The Big Rip singularity in (137) occurs since the potential is unbounded 
and goes to positive infinity when <f> — > —00. Eq.(139) tells that if V(4>) 
is bounded from above and has a maximum V m as for V{4>) = case, the 
energy density does not grow up infinitely and the Big Rip singularity does 



28 



not occur. We now assume for the large negative 0, the potential approaches 
a constant: 

V(<f>) — > V m (constant) when </> — > — oo . (144) 
In the region, Eq.(132) reduces to 

° = -@ + 3 "S< <"»> 

which can be solved as in (141). Then the energy density p^ (133) has the 
following form: 

P* = ^ + V ^- ( 146 ) 

The FRW equation becomes 

= 2^ + Vm ■ (147) 



K 



The first term in the r.h.s. could be neglected for a large universe. Then, for 
large a, one gets the deSitter space as a solution 

H ' 2 - h ■ (148) 

Thus, one way to avoid the singularity might be that, in the present universe, 
for large negative 0, there is an upper bound in the potential. 

Another way to argue is to take into account quantum effects, say, for 
conformally-invariant matter. Then the contributions coming from the con- 
formal anomaly to the energy density pa and pressure pa are (see [26]) 
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In general, with N scalar, N1/2 spinor, N\ vector fields, N 2 (= or 1) gravi- 
tons and TVhd higher derivative conformal scalars, b, b' and b" are given by 



b = 
b' = 



N + 6N 1/2 + 12iVi + 611iV 2 - 8N HD 
120(4tt) 2 

N + UN 1/2 + 62^ + 14117V 2 - 28iV H D 
360(4tt) 2 



b" = . (151) 



Near the Big Rip singularity, the scale factor a blows up, as in (137), 



at t — t s . Then the curvatures behave as R oc \t — t s \ and they become 
large. Since the quantum correction includes the square of the curvatures, 
the correction becomes large and important near the Big Rip singularity. 
Now the FRW equation has the following form 



6 



H 2 



7 
2 



dt 



+ V{<f>) + p A ■ 



(152) 



We now write H and <p as 

H = ho + Sh , 



^0 In 



t.-t 



(153) 



Here ho, t s , and £1 are constants. We assume that when £ — > t s , 5h and 
S(f> become very small compared with the first terms, respectively. In H, 
however, as the first term is a constant, only the second term contributes to 
and From (132) with the same exponential potential, one obtains 



dH 
dt 



= -7 

+o((t s -ty 1 ) 

which gives 

Vot\ = - 
Then from (152), we have 

37^000 



3h 



(t a -ty t a -tj o (£ s -£) 



+ 



2V tj 



l--8ct> 



1^1 



"2 



(154) 
(155) 







£.,-£ 



+ 6h (h + b"^j d ^ + ((t s -t)- 1 ) , (156) 
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and 



5h 



70o 



JX (t s ~ t) In 



t 8 -t 



2{lb + V) 

Here ti is a constant of integration. Since H — -, the scale factor is 



(157) 



a 





700 


i a -t 


4(§ 6+6 ») 


tl 





(158) 



In Sjgjr or there appears a logarithmic singularity. The behavior of the 
singularity, however, becomes rather milder than the case without quantum 
correction, where due to the singularity at t — t s , the universe cannot de- 
velop beyond the singularity. Since the singularity becomes mild due to the 
quantum correction and a and H are finite at t — t s , the universe can develop 
to the region t > t s . Then, essentially the Big Rip singularity is removed due 
to the quantum correction. 



5 Gravity assisted dark energy dominance 
and effective phantom/quintessence cos- 
mology. 

Recently, an effective phantom/quintessence description of the late time uni- 
verse was obtained via the introduction of a new higher-derivative coupling 
between matter and gravity [21]. It was shown that such a model may ex- 
plain the gravity assisted dark energy dominance. In the present section we 
will consider a simple (scalar) example of such kind of model when standard 
matter is also included. 
The starting action is: 

S = [ d 4 x^g{^R + R a L d + L m \ . (159) 



Here is the matter-like Lagrangian density (dark energy) and L m the 
Lagrangian density of the (standard) matter. By variation over g^ v , the 
equation of motion follows: 

= *? S =\ Hg^R - iff*") + + T% . (160) 
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Here the effective energy momentum tensor is defined by 

= -aR a - l W v L d + a (VV - g^V 2 ) + R a T^ . (161) 

and is given by 

V = ^Lj^ ( [ d 4 x^L d ) . (162) 



-g Sg, 



The standard matter part of the energy momentum tensor is also defined 

as 

^ s ^ii(J ■ (163) 

For simplicity, the Lagrangian density of a free massless scalar is considered 
as L d : 

L d = —d^d^ , (164) 

Note that for the above L d choice and with a higher derivative scalar curva- 
ture term in the gravitational sector (also without standard matter), model 
of this kind was discussed, with different purposes, in Refs.[22]. 
The equation of motion has the following form: 

= = -^Lfy (iTv^<TcU) . (165) 

v=g$(p v=g 

The metric is again chosen to describe a FRW universe with flat 3-space: 

ds 2 = -dt 2 + a(t) 2 ( dxt ) 2 ■ ( 166 ) 

i=l,2,3 

If one assumes that depends only on t ((f) = the solution of the 

scalar field equation (165) is given by 

= qa~ 3 R- a . (167) 

Here q is a constant of integration. Hence 

Q 2 

R a L d = — i— , (168) 
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which becomes dominant when R is small (large) compared with the Einstein 
term j^Riia > — 1 (a < — 1). Thus, one arrives at the remarkable possibility 
that dark energy grows to asymptotic dominance over the usual matter with 
decrease of the curvature. 

Combining (159) and (160), one gets S ~ f d^x^/^g j -^R + 2c f Ra j , 

_ 6 

which may indicate R ~ a q + x . Then the curvature R might be stabilized 
to have a non-trivial minimum due to the second term in (159). 

Substituting (167) into (160), the (/i, v) = (t,t) component of the equa- 
tion of motion has the following form: 



Pd 



3 

a 6 V / [ 4 4 

+ (\ + + a 2 ^j HE 2 + ^1 + ^ # 4 j . (169) 

Here p m is the energy density of the standard matter. Specifically, when 
a = -1, Eq. (169) looks like: 

If Pm = 0, this equation has only the trivial solution H = (a is a constant). 
When p m = 0, we can easily find the accelerating solution of (169) [21]: 

V 3t J • 3 (a + 1)" +1 (§ (2a - 1))° +2 

Eq. (171) tells that the universe accelerates, that is, a > if a > 2. Even for 
a < —1, by changing the time variable by t — > t — t (t is a constant), the 
universe is expanding and accelerating. In this case, however, there is a Big 
Rip singularity at t = t . 

For the matter satisfying the relation p = wp, where p is the pressure and 
p the energy density, from the usual FRW equation, one has 

aoct^^. (172) 



= ~^H 2 + p d + p, 
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For a oc t h ° it follows that 



+ (173) 



and an accelerating expansion (h > 1) of the universe occurs if 



-1 < w < ~ . (174) 



For the case of (171), one finds 

w = ^—^ . (175) 
1 + a 

Then if a < — 1, w < — 1, what corresponds to an effective phantom. In this 
case, changing t as to — t in (171), there appears a Big Rip singularity at 
t = t . In [23], however, it has been shown that the phantom energy with 
w < — 1 makes the radius of the wormhole spacetime (when it does occur) 
to increase in time and thus before the Big Rip the radius becomes infinite 
and, as a result, the Big Rip singularity may be avoided. 

It is interesting to investigate the stability of the solution in (171). For 
this purpose, we write the scale factor a as 

a = a t^ (1 + S) , (\6\ <1) . (176) 

Here a is given in (171). From (169), it follows that 

2 (2a - 1) (a + 1) (a - 1) 2 (2a 3 - 18a 2 - 33a - 7) d5 
*4 p ~dt 

9a(a + A)d 2 5 9ad 3 S 

The solutions of (177) are given as: 

5 oc t xo . (178) 

Here x is a constant. Then Eq. (177) gives 

= 9ax 3 + 9a (a + 1) + (-4a 3 + 27a 2 + 48a + 14) x 

-2 (a + 1) (2a - 1) (a - 1) . (179) 
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As clear from (175), when a — > +00, w > —1, which corresponds to 
quintessence and, when a — > —00, we find w < —1, which corresponds 
to a phantom. Imagine that parameter a — > ±00. If we assume x = 0(1), 
Eq. (179) reduces to 

= -4a 3 (1 + x ) + O (a 2 ) . (180) 

Hence 

x 1 . (181) 

With x = O(a), one obtains 

= -4a 3 x + 9a 2 x 2 + 9ax 3 + 0(a 3 ) . (182) 

As a result 

a 4a 

x ~0,-,-—. (183) 

Note that the first solution xq ~ corresponds to the solution in (181) since 
we have assumed x = O(a). The perturbation looks as 

5 ~ Srf- 1 + 5 2 t^ + 5 3 r^ . (184) 

Here #1,2,3 are constants. The second term in (184) may indicate the instabil- 
ity of the solution (171). When a — > +00, which corresponds to quintessence, 
the second term may become dominant when t — > 00. On the other hand, 
since the case a — > —00 corresponds to the phantom, we replace t — > to — t. 
The second term may become dominant near the Big Rip t — > to- This might, 
however, indicate that the Big Rip does never occur, since the Big Rip solu- 
tion is unstable. In other words, even if the present universe equation of state 
parameter looks as w < — 1, the universe might transit to another solution 
corresponding to the second term in (184). However, it is difficult to find the 
non-perturbative behavior of the solution corresponding to the second term 
in (184). 

The a — > — 1 case which corresponds to w — > — 00 can be also considered. 
Let us write 

a = -(l + e) , (e«l) . (185) 
It is natural to assume that e is positive. Eq. (179) gives 

= -9 (1 + e ) x 3 + 9ex 2 + (-3 + 18e) x + 12e + O (e 2 ) . (186) 
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Its approximate solution is 



x = 4e, -? c ±J=(l-^) . (187) 



The perturbation 5 is found to be given by 



S = 5 (t -t) 4e + (t - t)-'^ j<f c cos^= ^l-^ln (t - tfj 

+S s sm ( ( 1 - 7 -e ) ln(t„ -t) ) 1 . (188) 



Here 5o, S c , and 5 S are constants. The first term decreases near the Big Rip 
at t — to but the other terms oscillate rapidly near the Big Rip and the 
amplitude becomes large, which also shows the instability of the (transient 
acceleration) solution (171). This may be quite an acceptable result, with 
the condition that it lasts sufficiently enough to comply with observational 
data. The way to avoid a finite time future singularity due to the instability 
of the accelerating cosmology may deserve some attention. 
The case a > 1 corresponds to w < 0. With 

a = l + e, (0 < e< 1) , (189) 

Eq. (179) gives 

= 9 (1 + e ) x 3 + (18 + 27e) x 2 + (85 + 90e) x - 4e + O (e 2 ) . (190) 
Its solution is 

x = ^ + O(e'), -l±| + 0(e) . (191) 
Therefore 5 is given by 

S ~ 5 t§ + 1 j<5 c cos Q In tj + S s sin Q In tj | . (192) 

Here So, S c , and S s are constants, too. The first term increases with t, and 
thus the solution (171) may be unstable again. 

Even with numerical calculation for general a, there seems to be insta- 
bilities in the solution (171) if w < although the solution could be stable 
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if w > 0. This may suggest again that the acceleration of such dark energy 
universe might be transient. 

So far the standard matter contribution has been neglected. When p m ^ 
0, the simple assumption is that p m behaves as 

Pm = Poa' 13 . (193) 
Here p and (3 are constant. From Eq. (170) it follows 



t=xl^idaa^\l^ + %^. (194) 

For the case of a general a, if the dark energy density is dominant :pd ^> 
p m , the solution should behave as (171). Hence, the energy density p m of the 
standard matter evolves as 

_ o _ (a+l)/3 

p m = Poa t 3 . (195) 

On the other hand, pd behaves as 

Pd oc r 2 . (196) 

Since f3 — 4 for the radiation and (3 = 3 for the dust and the acceleration 
of the universe occurs when a > 2, we may assume (5 > 2 and a > 2. As a 
result, 

> 2 . (197) 

When time t grows, p m decreases further ias compared with p^. 

The essence of a gravity-assisted dark energy dominance is clearly seen in 
the example below. Let the standard matter dominate, as compared with the 
dark energy pm^ pd- Then the scale factor is given by that of the standard 
FRW equation: 

»=(^)*«»- m 

The energy density of the standard matter p m behaves as p m ~ r 2 . On the 
other hand, the dark energy behaves as 

p d ^ t 2a 'T . (199) 
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Hence if 

2a - — > -2 , (200) 

P 

the dark energy p d becomes larger as time passes. Eq. (200) can be satisfied 
if a, (3 > 2 as in the dark matter dominant case. 

Let us assume that in the early universe, the standard matter /radiation is 
dominant. The universe evolves according to (198). From (200), dark energy 
increases with time growth. When p d > p m , as in the present universe, 
the acceleration of the universe begins. Thus, in the future the accelerating 
universe evolves with the scale factor (171). However, that is most probably 
a transient acceleration. This is not strange owing to the fact that the above 
effective phantom/quintessence description is achieved by means of a higher 
derivative coupling between dark matter and gravity. It is known that higher 
derivative gravities (see [12] for a review) may have problems with unitarity 
(stability) when they are considered as fundamental theories. Hence, the 
model under discussion should be treated as a kind of effective matter-gravity 
theory. It would be interesting to analyze some astrophysical predictions 
of our model (say, rotation curves of galaxies) in the way it was recently 
discussed e.g. in [25]. 



6 Discussion. 

In summary, (phantom) scalar-tensor cosmology with an exponential scalar 
potential suggests the comely possibility of a dark energy universe with an 
equation of state parameter w which is negative and very close to —1. Con- 
venient choices of the theoretical parameters may shift the explicit value of w 
from above to below —1. Moreover, such a universe naturally admits a (tran- 
sient or eternal) acceleration phase. A very nice property of this theory is that 
what appears as a phantom in one frame may look as a standard scalar in an- 
other frame. It is demonstrated that in the situation when a finite time future 
singularity is predicted by the growing phantom energy density, the consid- 
eration of quantum gravity effects might drastically change the future of our 
universe, removing the singularity in a quite natural way. From another side, 
it is also shown that a higher-derivative gravity-matter coupling term being 
not the phantom may in fact provide an effective phantom/quintessence de- 
scription of the late-time universe, suggesting the possibility of a dark energy 
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model of a brand new type. In this case, gravity makes dark energy to become 
the (evolving) main contribution to the total energy density — as compared 
with the standard matter/radiation, which was initially dominant — what 
leads to the apearence of the phase of transient acceleration. 

Current attention to phantom models as dark energy candidates is not 
driven by the internal consistency and/or beauty of this theory, which still 
contains a number of non-fully-resolved problems, as we have already men- 
tioned. Rather, it is the lack of a good theoretical understanding of the 
present universe coming from more usual theories what calls for alternative 
explanations to be considered, on one hand. From another, one sees also that 
the cosmic zoo structure which emerges from these alternatives is so rich and 
suitable at times, that some concepts which so far seemed to be strange (like 
the idea of negative energy itself) deserve to be investigated with care and to 
the end. In this respect, even the mild indications which have been reported 
of a possible phantom origin coming from string/M-theory or on the chance 
(described above) to avoid the Big Rip catastrophe by taking properly into 
account the quantum effects seem indeed very promising. 

Acknowledgments 

We thank James Lidsey, Alexei Starobinsky and Paul Townsend for fruitful 
discussions. This investigation has been supported in part by the Ministry 
of Education, Science, Sports and Culture of Japan under grant n. 13135208 
(S.N.), by the RFBR grant 03-01-00105 and LRSS grant 1252.2003.2 
(S.D.O.), and by SEEU and DGICYT (Spain), grant PR2004-0126 and 
project BFM2003-00620, respectively (E.E.). 

A Appendix 

In this Appendix several remarks are made about the possible origin of the 
phantom-related models coming from the higher dimensional theories con- 
sidered in the paper. This topic was widely investigated in the Kaluza-Klein 
context (in relation with the string/M-theory), and for a recent discussion 
the reader is addressed to [17]. We start from a 4 + n-dimensional spacetime, 
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whose metric is given by 



n 



ds 2 = 9^dx^dx u + e 2 ^ M) ^ il,;A'< J ■ 



(201) 




For simplicity, one may assume that the metric gij corresponds to an Einstein 
manifold, where the Ricci tensor Rij constructed from g^ is proportional to 
gif = kc/ij. Here A; is a constant. When n — 1, k always vanishes 
(k = 0). When n > 3, the above metric is given as the solution of the n- 
dimensional Euclidean Einstein equation. When n = 2, since the 2d Einstein 
equation is trivial, in the conformal gauge the above condition for the Ricci 
tensor is the Liouville equation. Under the above assumptions, the 4 + 
n dimensional Einstein action with matter field x (bosonic sector of some 
higher-dimensional supergravity) can be written as 



Here V n is the volume of the n-dimensional manifold whose metric tensor is 
given by cjij. We should note that the kinetic energy of 4> becomes negative, as 
for the phantom. Rescaling the 4-dimensional metric g^ by g^ v — > e - ™^^^, 
the action (202) can be rewritten as 




(202) 




(203) 



Now the kinetic energy of cf) is positive. If we further rescale <f) by 




(204) 
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the four-dimensional action looks like 



(205) 



S 4+n = ^ J ' d A x^-g(R-^d^d uV 

+nke-^^ - U^X ~ e^^U( X )^ . 
The above action belongs to the same class as (43) by identifying 

W(^x) = e~ ip ^U( X ) ■ (206) 
Comparing the above expression with (47), one obtains 



3n 2 + rj 



n + 2 ip 
Then, if 



(207) 



/n + 2 

0<^o<2^/^ r , (208) 

Eq. (53) follows. If 

U(x) = W e^ x , (209) 

we can obtain for w a value less than —1, and correpsondigly the universe 
expands with acceleration, by a proper choice of the parameters. 

One may consider the product compactification to be more general than 
(201): 

ll,v=0 i,j=l 
N 

+e ^n J2 g?}dCd? . (210) 
I,J=1 

Since the scalar curvature R( d + n+N ) in (d + n + iV)-dimensional spacetime is 
given by 

R (d+n+N) = R (d) + k (D e -2<pW _ + ^^0(1)^0(1) _ 2nV 2 (1) 

+jb (2) e -2* (2) _ N (N + 1)^ 2 W 2) _ 2iVV 2 (2) 
-2niVa M (/) (1) a At ( /) (2) , (211) 
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if we start from the (d + n + iV)-dimensional action coupled with the scalar 
field x with potential U (x) , we obtain 

1 



S (d+n+N) = J_ f d d+n+N x /_ (d+n+N) 
K 2 J 



R {d+n+M) -^d, X d»X-U{x) 



v n v N 



-iVa M (/) (2) ^0 {2) + d M (n0 (1) + iV0 (2) ) (n0 {1) + iV0 (2) ) 
+A ,(i) e -2^( 1 ) + fc (2) e -20p) _ \ d ^ x _ c/( x )l . (212) 

If is a constant, one may regard it as the cosmological constant. Further 
rescaling the metric tensor by 



9/w 



9»v 



(213) 



we get 

g(d+n+N) 



V n Vj 



N 



K 



d-2 



d d xy/-g {d) R (d) - nd^dr^ - Nd^d^^ 



+k (%-M n <t> (1) + N <t> i2) )- 2 <t> (1) + A .(2) e -^K 1 )+^< 2 ))-2^) 
1 



(214) 



In this frame, the kinetic term of the matter field x does not directly couple 
to the scalar fields <fi^ and (p^ . Then the Newton law is not violated in the 
leading order of perturbation. The kinetic terms of the fields (fP and (p^ 
can be diagonalized by 



(jP) = y-(j)+ + y+ < fT , 0( 2 ) = 



Here 



1 



1 ± 



(d-2) 2 (n-N) 2 (l + ^y 
An 2 N 2 



(215) 



(216) 
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Hence 



+ A;( 1 ) e -{(^+ 2 )^ + ^ + }^ + -{(^+ 2 )^-S^}^ 



-\q )xX q» x _ e ~ih{( n y-+ N y + )t + +( n y + - N y-)<t>-}u(x)\ (217) 



Here, 



x 



± 



If Ar n 2 + N 2 r-\ 



Ar n 2 + N 2 \ 2 ( n + N 

D = [n + N + — - AnN 1 + 



d-2 / V d ~ 2 

,^2 f n + N\ 2 4n 2 N 2 . , 

= (n — TV) ( 1 H — : — - + ^>0. 218 

V ; V d-2 J (d-2) 2 K J 

For simplicity, the case = U(x) = X — and d — 4 is considered. 
Comparing the action (217) with (43), we may identify 



'2x+ 

(p <-> y — — + , /tx^V2a;"0~. (219) 



From (47), one sees that 

Wo = -A; (1) 

2 + 7] 



{(n + 2)y- + Ny + }^: 



tp L J V 2x + 

(n + 2)y + - Ny- 
Xo V2x~ 

Using (50), an expression for 77 can be found: 



(220) 



= 2_ ^ o;-{(n + 2) j/ - + iV ?/ +} 2 

1 + C' x+{(n + 2) ?/ +-iV ?/ -} 2 ' 
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Since x ± > (218), £ > and therefore rj satisfies (53). Then w > — 1 and 
there is no phantom. More explicitly w is given by 



, 2 {(n + 2)y- + Ny+] , nnn . 
w = -l + — r JLJ ~ • 222 

Numerically, with n = 2 and N = 5, w = 0.00269554; with n = 3 and iV = 4, 
id = —0.0586691. Although not realistic, if we choose n = 6 and iV = 31, it 
follows that w = -0.353435. 
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